Optimizing the second hyperpolarizability with minimally-parametrized potentials 
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The dimensionless zero-frequency intrinsic second hyperpolarizability jim = 7 / AE±q 'm~ 2 (eh) 4 
was optimized for a single electron in a ID well by adjusting the shape of the potential. Optimized 
potentials were found to have hyperpolarizabilities in the range —0.15 J$ ji n t J$ 0.60; potentials 
optimizing gamma were arbitrarily close to the lower bound and were within ~ 0.5% of the upper 
bound. All optimal potentials posses parity symmetry. Analysis of the Hessian of 7i n t around 
the maximum reveals that effectively only a single parameter, one of those chosen in the piecewise 
linear representation adopted, is important to obtaining an extremum. Prospects for designing new 
chromophores based on the design principle here elucidated are discussed. 



I. INTRODUCTION 

Developing materials with high electronic nonlinear 
susceptibilities is of fundamental importance for a wide 
variety of applications such as optical solitons, phase con- 
jugate mirrors and optical self- modulation [IJ [2] • These 
susceptibilities are defined by considering a material in 
the presence of an electric field E and expanding the in- 
duced polarization in a Maclaurin series, 



P = aE + PEE + jEEE + 0(E 4 ), 



(1) 



where the susceptibilities a, f3 and 7 are generally 
frequency-dependent tensor quantities. Here, we con- 
sider the zero-frequency limit of these quantities in the 
non-resonant regime i..e. where all frequencies are much 
less than any resonant frequency of the electrons and mo- 
tion of the nuclei is neglected. A remarkable result due 
to Kuzyk [5] is that quantum mechanics requires that 
the first and second hyperpolarizabilities (3 and 7 are 
bounded: specifically, 7 obeys the inequality, 
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where N is the number of electrons, E\$ is the energy 
difference between the ground and first excited states and 
in is the electron mass. It is natural to define the intrinsic 
hyperpolarizability as a figure of merit to characterize the 
proximity of a given system to this limit, 



Itnt = 7/7o" 



(3) 



and to ask: how to create materials that achieve opti- 
mal 7i„t? The discovery of the bounds ^ has moti- 
vated a number of experimental studies that have demon- 
strated that carefully tuning the electronic states and 



geometry of chromophores can lead to higher second 
hyperpolarizabilities jl]-[7]. Generic design principles mo- 
tivated by fundamental theory would therefore be desir- 
able. Unfortunately, the procedure used to derive the 
bounds Q cannot directly provide these; they were ob- 
tained by optimizing 7 for a three-level ansatz with re- 
spect to the dipole matrix elements and energy level spac- 
ings E — E10/E20 and not by constructing an explicit 
potential. Indeed, the assumptions behind the deriva- 
tion have been questioned[8, 9 and the limits need not 
be achievable with a local potential; it has been specu- 
lated recently these may require exotic HamiltoniansfTU] . 

Subsequent work, following approaches developed in 
earlier studies of B\ll\ IT?] , has attempted to address 
this in two ways: First, by identifying universal features 
of Hamiltonians near the fundamental limit by a Monte 
Carlo search[T5] and, secondly, by numerically optimiz- 
ing with respect to the shape of a local potential) 10|. 
This latter work found potentials which have second hy- 
perpolarizabilities in the range —0.15 < ji n t < 0.60, 
which represents an apparent bound that is more re- 
strictive than the bound of Moreover it was demon- 
strated that the optimized potentials spectra and dipole 
moments were broadly consistent with those identified in 
the earlier Monte Carlo study. 

While these strategies provide useful goals for chemists 
attempting to design new nonlinear chromophores, they 
do not provide insight into which features of the potential 
are necessary to optimize 7j n4 , or how many free parame- 
ters should be necessary to achieve optimal or near opti- 
mal 7. In a previous paperfll]. we developed a technique 
to examine the analogous question for (3: by optimizing 
potentials described by increasing numbers of free param- 
eters and examining the eigenvalues of the Hessian matrix 
at each maximum, we identified the combinations of pa- 
rameters most important to the optimization. The anal- 
ysis revealed that effectively only two parameters were 
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necessary to maximize /3, and hence that a surprisingly 
broad range of potentials with high j3 exists around each 
maximum. 

In this work, we apply the same technique to the prob- 
lem of optimizing 7. At first sight, the problem appears 
to be more difficult than that for j3 since the expression 
for 7 is much more complicated and the bounds for pos- 
itive and negative 7 are different. Remarkably, however, 
we will show that effectively only one parameter is neces- 
sary to optimize 7 in either direction and, moreover, that 
in each case it is one of the parameters utilized in our rep- 
resentation of the potential. In this sense, we are able to 
suggest much more clearly a possible design strategy for 
materials with high 7 than for /3 within the limitations 
of the model. At least within our representation of the 
potentials, we find that the potential that maximizes 7 
is rapidly varying, while for negative 7, as for /3, quite 
generic, slowly varying potentials are adequate. The pa- 
per is organized as follows: in section II the calculations 
performed are described; the results are presented and 
discussed in section III; conclusions are drawn in section 
IV. 



II. MODEL 

It is first necessary to generalize the method described 
our previous paper on optimizing the intrinsic first hy- 
perpolarizability /3i n t |14j : in the present work, 7^ is to 
be optimized by adjusting the shape of a one-dimensional 
piecewise-linear potential. Such a potential with N + 1 
segments may be represented, 

{A x + B Q x < x 
A n x + B n x n -i < x < x n , n € {1, ...N - 1} 
A N x + B N x>xn-Xi 

(4) 

with the positions x n and slopes A n as the adjustable pa- 
rameters and where the B n are chosen to enforce conti- 
nuity. Because 7^ is invariant under trivial translations 
and rescalings of the potential, some of these parameters 
were fixed xq — 0, Bo = B\ = 0, and A\ = ±1. These 
choices, together with a change of origin and rescaling 
allow for any potential. Thus maximizing with these 
constraints allows faster optimization. Furthermore, the 
left- and right-most slopes are required to be negative 
and positive, respectively, ensuring only bound electron 
states. Finally, for technical reasons, having to do with 
the asymptotic behavior of the Airy functions introduced 
in eqn[7]below, it is difficult to allow the sign of any slope 
to change during an optimization. In consequence, we 
have chosen to restrict | A4 |> .005, and, as appropriate 
to do separate optimizations for each interesting sign of 
each slope. 

A second representation for the potential was also con- 
sidered where parity symmetry was specifically enforced. 
This was motivated by previous work [15. which identifies 
parity as important for optimizing 7mt, particularly for 



the lower bound. The potentials with enforced V sym- 
metry were constructed on the half line x > with N 
segments, 



V(x) = 



A n x + B Ti 



i_l < x < x n , n e {1, ..JV - 1} 



Ajqx + Bn x > x^-i 



( 5 ) 

with xo — and requiring V(—x) — V(x). Again x n and 
A n are adjustable parameters and xo = 0, Bo = B\ = 
are fixed. The parameter A\ was set to either —1 or +1 
to study the consequences of both cases. 

For such a potential with a uniform applied elec- 
tric field of strength e, the wavefunction obeys the 
Schrodinger equation in each segment, 



l_d*_ 
~2dx 2 



+ (A n + e)x + B ri 



1p n = Elp n , 



(6) 



in units such that e = 1, K = 1, and m e = 1. The 
solution in each segment is written in terms of the well- 
known Airy functions, 



C n Ai 



^2(B n - E + x{A n + e)) 
{A n + e) 2 /3 

V2{B n - E + x(A n + e)) 
(A n + e)2/3 



(7) 



To solve for the coefficients C n and D n the usual bound- 
ary conditions are imposed, i.e. that the wavefunc- 
tion ip(x) and its derivative ip'( x ) are continuous at the 
boundary between segments. Additionally, in the end 
segments, the wavefunction must vanish as x goes to ±00 
fixing Dn = 0. There are a total of 2N linear equations 
in the coefficients for the arbitrary case and 4iV — 2 equa- 
tions and coefficients for the "P-symmetric case, which can 
be written in matrix form 



W-u = 



(8) 



where u is a vector comprised of the C n and D n coeffi- 
cients and IF is a matrix which depends on £, £ and the 
parameters A n and x n . 

The allowed energy levels are found by numerically 
finding the roots of 



det W = 



(9) 



with e = 0. It is readily possible, as previously done 
for /?, to obtain from Q an expression for the second 
hyperpolarizability, 
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this is achieved by repeatedly differentiating the matrix 
W using the Jacobi formula, 
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det IF = Tr 



(11) 
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where adjVF is the adjugate of W (since W is singular), Having performed similar calculations to those in [14], 
and applying the chain rule we arrive at the expression 

dW _ dW dWdE 

~d7 ~ Ih + Imlk' ( ' 



d 4 E 
~d^ 
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(13) 



where 



W"' = 



d 4 W dW d 4 E 



de 4 dE de 4 



(14) 



From this, 7 is readily obtained and the intrinsic second 
hyperpolarizability ^ int — 7/7™ ax calculated for a given 
set of parameters. 

The quantity ji n t was optimized numerically for both 
arbitrary and "P-symmetric potentials with varying num- 
bers of segments using the FindMaximum function of 
Mathematica 8, an implementation of the Interior Point 
method for constrained optimization. Both maxima and 
minima of 7;„ t were obtained from a large number of 
randomly generated starting points, and also manually 
chosen starting points with large values for 7i n j. Once 
an optimum 7i„ t was found, the extent to which each 
of the parameters was important to the extremum was 
characterized by calculating the Hessian matrix, 
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where the Pj are the parameters, and calculating its 
eigenvalues and eigenvectors. Since the Hessian matrix 
characterizes the local curvature of the objective func- 
tion in the parameter space around the extremum, these 
quantities give the magnitudes and directions of the prin- 
cipal curvatures. As stressed in previous work[14 these 
curvatures implicitly depend on a measure implied by 
this equation that is peculiar to our numerical parame- 
terization of the problem. More physically relevant mea- 
sures can also be used to calculate eigenvalues. While 
these make quantitative changes in the eigenvalues and 
vectors, they do not make qualitative changes, and we 
give results for this "numerically natural" measure below. 



III. RESULTS AND DISCUSSION 

Our optimized potentials, together with the ground 
and first excited state wavefunctions, are displayed in 
Fig. [I] for both arbitrary (eq. [4]) and "P-symmetric (eq. 
[5]) parametrizations. The associated parameter values 
are listed in Table [I] The potentials and wavefunctions 



are displayed on a transformed position and energy scale, 

x = (x-ix^/iEi-Eo) 1 / 2 
V(x,{P}) = (V(x,{P})-E )/(E 1 -E ) (15) 

such that the ground state energy is Eq = 0, the differ- 
ence between the ground and first excited state energy is 
E\ — Eq = 1 and the position expectation value for the 
ground state is (xq) = 0. This rescaling does not change 
"fint and permits convenient comparison of the results of 
each optimization. To identify the relative importance of 
each of the parameters to the optimization, the results of 
the eigenanalysis of the Hessian matrix are also displayed 
for selected potentials in Fig. [T] the j-th eigenvalue of 
the Hessian, hP, is listed alongside a plot of the variation 
in the potential AV^x) in the direction of the associated 
eigenvector, 



AV j (x) 



dV{x,{Pi+a^}) 



(16) 



where vj is the i component of the j-th eigenvector. 
Note that the values of V and x in the right hand side 
of (16 1 are renormalized as a function of a using (15 1 so 



that the variations presented automatically preserve the 
properties (xo) = and E\ — Eq = 1. 

The first set of results displayed in Fig. [l|a) are op- 
timized potentials with no enforced symmetry and spec- 
ified by 3 or 5 free parameters. The optimized ji n t for 
both the lower and upper bounds of ji nt potentials fall 
somewhat below the apparent bounds observed in |13| . 
For the upper bound, the best results assuming all slopes 
except the first are positive are what appears to be a lo- 
cal maximum value of 7, nt ±2 0.43 [Fig. [lja)(ii)]. A 
similar value was also found in [10] for potentials with 
no constrained symmetry, but the potentials found here 
do not closely resemble those found in that work. The 
eigenvalues and eigenvectors displayed below the poten- 
tial in [Fig. [l|a)(ii)] show that only two of the eigenvalues 
are significant in magnitude and are associated with the 
shape of the potential in the middle while the small eigen- 
values are associated with the outer slopes. These results 
are reminiscent of those found for the first hyperpolariz- 
ability, where /3i„ t was found to approach its maximum 
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Bound 


Description 


Ket. 


-fint 


A A 


A 


X\ 


X2 




3 param. arb. 


(a)(i) 


0.43277 


2.33932 0.10189 





1.04960 





Upper 


5 param. arb. 


(a)(ii) 


0.43602 


4.28494 0.08829 


0.37688 


1.42131 


3.95728 


2 param. V 


(b)(i) 


0.58220 


0.00550 




0.57426 






3 param. V 


(b)(ii) 


0.59707 


0.00500 




0.88533 


7.95755 


Lower 


3 param. arb. 


(a)(iii) 


-0.11409 


82.471 237.57 




1.25227 




2 param. V 


(b)(iii) 


-0.14996 


31.6157 




1.48005 





Table I: List of parameters of optimized potentials 
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Figure 1: Optimized potentials for the upper bound of 7i„t 
with (a) no enforced symmetry and (b) P-symmetry. The 
energies of the ground and first excited state are indicated by 
horizontal lines; the corresponding wavefunctions are also dis- 
played. The plots have been rescaled to facilitate comparison 
by ensuring E 1 — Eq — 1 and (xq) — while preserving 7i„t. 
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Figure 2: Optimized potentials for the lower bound of 7i„t 
with (a) no enforced symmetry and (b) P-symmetry. 
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Figure 3: Five parameter potentials with no enforced symme- 
try optimized using two-parameter 'P-symmetric optima as 
the starting point, (a) Upper and (b) lower bounds. 



value for the same class of potentials with a similarly 
small number of parameters and analysis of the Hessian 
revealed that only effectively two parameters were im- 
portant to the maximization. For the lower bound, the 
3 parameter system [Fig. |2^a)] converges on a shape ap- 



proaching a square well. The 5 parameter system also 
converges on a P-symmetric potential and, because of 
this, further discussion of the lower bound is deferred to 
a subsequent paragraph. 

A possible explanation for the fact that the optimized 
positive ji n t for arbitrary potentials in Fig. [T] falls short 
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of the bounds established in [3J [TOJ H3| is that the hyper- 
polarizability is sensitive only to some features of the 
potential and that many local extrema exist. Rather 
than making extensive runs starting from a variety of 
potentials, we chose to use the more constrained subset 
of ^-symmetric potentials (|5| and found that, indeed, 
even with only 2-3 parameters, much higher values of 
7i„i could be obtained as shown in Fig. [ljb). This is 
reminiscent of the observation in [TU] that these bounds 
could only be reached if a ^-symmetric starting point 
was used for the search; in this work, we not only enforce 
the symmetry of the starting point but at all times in the 
optimization. 

For the upper bound, we found a two-parameter po- 
tential close to the apparent maximum but below the 
theoretical maximum [Fig. [T[b) (i)]. The shape is char- 
acterized by shallow outer slopes with a divot in the cen- 
ter; the ground state wave function is localized to the 
divot with the highly delocalized first excited state above 
the divot. It was found to be necessary to constrain the 
slope A2 > 0.005 since the method of calculating ji nt pre- 
sented in section [TT] fails for shallow slopes. To avoid the 
unphysical feature of delocalized higher excited states, a 
distant wall was added to construct a three-parameter 
potential [Fig. [ljb)(ii)]. Such a change is expected a 
priori from previous work|16j to make no significant dif- 
ference to ji n t as it is far from the region where ijio and 
ipi are large. For this well-like potential, we performed 
a maximization, adjusting A2, x\, and X2 while fixing 
the outer walls to have a large slope (A3 — 100). The 
best potential found possesses ji nt = 0.59707, and has 
A2 = 0.005 which is the shallowest slope allowed by the 
constraint. The eigenvectors of the Hessian, calculated 
for the subspace excluding A 2 , are well aligned with the 
parameters: the most significant eigenvector corresponds 
to Xx, the outer boundary of the divot. The other eigen- 
vector corresponds to X2, the position of the outer walls; 
this would be expected to have relatively little influence 
on 7i„ t since it controls a feature where the ground state 
wavefunction is small. Fixing the outer slopes at A3 = 100 
as above, we also attempted to optimize ji nt with A 2 con- 
strained to be negative. It was found that ji nt increased 
as the slope approached zero, until a point was reached 
where the calculations became numerically unstable due 
to the asymptotic properties of the Airy functions. The 
highest value of 7s n j which was found within a region of 
parameter space for which the calculation was still stable 
was 7j„i = 0.5915, lower than the current maximum. We 
then performed similar optimizations on potentials with 
strict hard wall boundary conditions. The calculations 
for these potentials were numerically stable in all regions 
of parameter space which were explored: For A2 > 0, 
a maximum of "f int = 0.5959 was found. For A 2 < a 
higher maximum of 7^ = 0.5968 was found, though this 
is still lower than the current maximum of 7^ = 0.5971. 
Since a higher 7^ is found in potentials with A2 < 
than for potentials with positive A 2 in cases with strict 
hard wall boundary conditions, we speculate that a value 



of 7mt higher than the current maximum found might be 
found for the finite A3 case. Nonetheless, we do not ex- 
pect to see a significant improvement as the current max- 
imum is already within ~ 0.5% of the maximum value 
found in previous studies. 

For the lower bound, a potential with the best value 
of lint = —0.1500 was found using only two parame- 
ters [Fig. [2^b)]. This potential is characterized by steep 
outer walls and a "bump" in the middle: the ground state 
and first excited state wave functions cover the same spa- 
tial extent, but the bump causes the ground state to be- 
come spread out and relatively flat. Eigenanalysis of the 
Hessian of 7^ about this solution shows that one eigen- 
value is significantly larger than the other, indicating that 
only one of the parameters is physically relevant. More- 
over, the eigenvectors of the Hessian for this potential 
are aligned with the parameter space chosen to represent 
the potential. The higher eigenvalue is associated with 
an eigenvector along the x\ direction, which determines 
the position of the outer walls; the smaller eigenvalue is 
associated with the parameter that controls the slope of 
the outer walls. Because j int is invariant under rescalings 
of the form (15 1, a potential with identical 7i„t can be 



constructed for a well of arbitrary width by tuning the 
slope of the bump. 

The two parameter "P-symmetric potentials identified 
as the satisfying the apparent lower bound can be equiva- 
lently represented by a five parameter arbitrary potential. 
Optimization of the five parameter potential indeed finds 
this potential as the apparent global maximum. There 
are therefore no directions in this new parameter space 
which lead to a higher 7i„ t , despite relaxing the require- 
ment that x_„ = —x n and A_„ = —A n . While the 
existence of asymmetric potentials with more negative 
7i„t cannot be ruled out, our analysis confirms that a V- 
symmetric potential satisfies the apparent lower bound. 
We repeated this procedure for the upper bound using 
the two parameter 'P-symmctric potential displayed in 
Fig. [2jb) as the starting point for optimization. It was 
found that despite relaxing the symmetry constraint no 
further improvement could be made so that the potential 
of Fig. [2^b) is also a local optimum with respect to the 
expanded parameter space. Analysis of the hessian was 
also performed on both of these five parameter optima. 
The eigenvalues and associated eigenvectors, shown in 
Fig. [3j reveal that for both upper and lower bounds 
the variation in the potential associated with the largest 
eigenvalue is indeed asymmetric. 

The fact that we are able to achieve optimized j int 
within 1% of previous limits with far fewer parameters 
than in other representations of the potential[l(j], and 
also fewer parameters than required to numerically opti- 
mize Pi n t, is surprising since the calculated expressions 
for ji n t are far more complicated than those for /3j n {. 
Since the dimensionality of this parameter space is so 
small, it is instructive to visualize it directly: For our 
two-parameter optimizations, the value of 7^4 is plot- 
ted over a portion of the parameter space [Fig. |4ja) and 
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Figure 4: Visualization of the variation of ji n t as a function 
of the parameters xi and A2 of a two-parameter P-symmetric 
potential for the (a) upper bound and (b) lower bound. The 
region of parameter space where 7i nt is within 2% of the max- 
imum value is highlighted in (a) and indicated by an arrow. 
jint is also shown in the (E, X) parameter space for the (c) 
upper bound and (d) lower bound. 



(b)]. The plot for the lower bound merely illustrates the 
results obtained from analysis of the Hessian, i.e. that 
the minimum is strongly curved about the optimal value 
of x\ but shallow with respect to A^. The plot for the 
upper bound is more interesting: the region of parameter 
space for which ji nt is within 2% of the maximum value 
obtained is highlighted showing a clear ridge. 

By calculating values of X = xqi/xq^ and E = 
E01/E02, the natural parameters of the three-level 
ansatzpj], as a function of (A2,x\), we are able to dis- 
play 7 ijlt re-parametrized in (E, X) space [Fig. |4jc) and 
(d)[. Here, x^\ ax = l/y/2E^ in our units. Notice that 
the entire region explored in the numerical parametriza- 
tion (A2,xi) collapses onto a narrow, elongated region 
in (E, X) space for the upper bound [Fig. mc)] and a 
complicated curved line for the lower bound[Fig. [4jd)[. 
These plots confirm the results of eigenanalysis of the 
Hessian: that essentially only a single parameter (or 2 
for the upper bound) characterizes optimal 7m*. 

The results of the three-parameter optimization and 
the plot in fig. Q both suggests that the truly optimal 
■p-symmetric potential for the upper bound has shallow 
outer slopes A2 — > 0. Such a potential can be trans- 
formed, using ( |15| |, to a potential of equivalent 7i„ t but 
where the outer slope is unity and the central well is far 
narrower and sharper. Since the central divot for the 
transformed potential resembles a Dirac delta function, 
we studied the second hyperpolarizability of the family 
of potentials 
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Figure 5: Plot of j int , E = E 10 /E 2 o and X = x 10 /x™? x as a 
function of a for the for the potential V = \x\ — aS(x). 



when a = 1.69552. Despite the simplicity of this one- 
parameter potential, the result is only 3% smaller than 
the best reported so far and only fractionally smaller than 
those found with the two and three parameter potentials 
above. 

We now turn to the question of whether the simple, but 
optimal or nearly so, potentials that were identified above 
resemble either those previously found [10 or possess the 
universal features identified in [13]. Comparing our best 
optimized potentials and wavefunctions to those in |10| . 
some qualitative similarities are apparent. For the upper 
bound, the potentials in Watkins et al. are roughly sym- 
metric near their lowest point. They feature a central 
divot within a wider well where the ground state wave 
function is localized within the central divot and the first 
excited wave function is relatively delocalized compared 
to the ground state. For the lower bound, the potentials 
feature a steep well within which both the ground state 
and first excited state are localized, and a central bump 
which causes the ground state to be spread out within the 
well. These qualitative features are shared by the poten- 
tials obtained in this work, but many extraneous details 
are removed by the highly constrained, judiciously chosen 
representation. 

In table|IT]we display values of X and E, together with 
the values of these parameters that extremize 7^4 for the 
three-level ansatz|3] and those possessed by the best pre- 
viously found potentials[10j. Values are also shown for 
two elementary potentials, the triangle well and infinite 
square well. The results for the upper bound are quite 
consistent with those of Watkins et al. if the breadth in 
the range of these values identified by the Monte Carlo 
study [T3"] is taken into account. The values of E and 
X are also displayed in Fig. [5] as a function of a, the 
strength of the 6 function, in the potential ( 17 1; E is a 



where a is the single adjustable parameter. Values of 7^ 
as a function of a are displayed in fig. [5] and the max- 
imum value is found to be 7j ni = 0.58194 which occurs 



monotonically increasing function of a while X is a mono- 
tonically decreasing function. Crudely, these explain the 
existence of a maximum 7i„ t as representing the trade-off 
between increasing the motion of the electron (associated 
with high X) versus enhancing transitions to other states 
(associated with low E). 
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Bound 



Potential 



Jint 



E 



X 



Upper 



3-level ansatz[3 
Best from [TUj 
Arb. 5 param. 
V 3 param. 
\x\ — a8 
Triangle well \x\ 



1 

0.5993 
0.43601 
0.59707 
0.58194 

0.1392 





0.5767 
0.5765 
0.8665 
0.6842 
0.5918 





0.5150 
0.6512 
0.7468 
0.7385 
0.9911 



Lower 



3-level ansatzg] =025 ±1 
Best from[T0] -0.1500 0.1493 0.6658 
V 2 param. -0.1500 0.2855 0.9416 
\x\ - a8 -0.1236 0.2438 0.9222 

Infinite square well -0.1262 0.3750 0.9801 



Table II: Second hyperpolarizabilities and physical parame- 
ters X = xoi/xqi 1 * and E = i?oi/-Eo2 for the optimized po- 
tentials obtained in this work. 



IV. CONCLUSION 

We have optimized the intrinsic second hyperpo- 
larizability j int of a piecewise linear potential well 
with respect to parameters that control the shape of 
the potential. We found solutions that lie within 
the range —0.15 < ji n t J$ 0.60 in agreement with 
the apparent bounds established in previous numerical 
optimizations^!)]; these both fall short of the Kuzyk 
limits[5]. By using two types of potential, one where 
all slopes were allowed to vary and another with ex- 
plicitly enforced symmetry, we demonstrated that V- 
symmetric potentials satisfy the apparent lower bound 
for 7i„ t and come within ~ 0.5% of the apparent upper 
bound. The parametrization used constrains the poten- 
tial to be smooth, preventing the occurrence of rapid 
oscillations which do not affect 7^ [IBJ. Because of 
this and the strong symmetry constraint, the optimal V- 
symmetric potentials found were characterized by only 
2 — 3 parameters. Of these, a posteriori analysis of the 
Hessian revealed that effectively only one or two, for the 
lower and upper bound respectively, were important 

These results are reminiscent of those obtained ear- 
lier for /3i ri t[14J, yet the number of parameters required 
to optimize 7i, it appears to be smaller even though it is 
a more complex object, containing more terms and in- 
volving higher derivatives. At least part of the reason 
for this is that for ji, lt there exists a "compatible" sym- 
metry operator, V, which can be used to constrain the 



shape of the potentials; this was not possible for flint 
where 'P-symmetric potentials automatically have fl = 0. 
However, even though we have shown that "P-symmetric 
potentials can have optimal or near-optimal second hy- 
perpolarizabilities, it is not clear whether the apparent 
upper bound 7^ ~ 0.6 achieved by Watkins and Kuzyk 
can be achieved with a "P-symmetric potential or whether 
a small amount of asymmetry is necessary. Moreover, 
the reason why local potentials fall short of the Kuzyk 
bounds remains opaque. 

The small parameter space allows us to propose a clear 
design paradigm for new chromophores, within the limi- 
tations of the model one-electron ID system studied. Ne- 
glected here, for example, arc multi-electron interactions, 
molecular ordering and inter-molecular electron hopping. 
Nonetheless, the potentials obtained could be realized, 
for example, by a centrosymmetric molecule with a cen- 
tral attractive or repulsive group — for positive or nega- 
tive 7j„ t respectively. The strength/electronegativity of 
the central group and the ratio of the length of the cen- 
tral and peripheral groups can then be tuned to give high 
"fint- Unfortunately, most practical chromophores are 7r 
conjugated systems in which there are approximately as 
many electrons as there are sites on the molecule, thus 
the consequences for them from this single electron cal- 
culation are clearly very speculative. Nevertheless, ho- 
mologous sequences already studied for high 7 as a func- 
tion of chain length, e.g. [BJ could likely be enhanced by 
including such a central group with different electroneg- 
ativity. Any other means of achieving a potential well or 
inducing a significant phase shift in wavefunctions pass- 
ing through the center of the molecule is also likely, even 
in multi-electron systems, to offer a route for achieving 
larger "fi n t- The present analysis provides other impor- 
tant insights: first, that since the "true" parameter space 
for -fi n t is so small, only rough tuning of the molecular 
design ought to be necessary. Secondly, this work again 
confirms that there are a large set of modifications to op- 
timized potentials, e.g. rapid oscillations, that will not 
change 7^ and need not be considered in planning what 
molecules to synthesize. As has been previously noted, 
materials with high j int could also be realized more di- 
rectly in other ways, such as through composite materi- 
als. 



[1] R. W. Boyd, Nonlinear Optics, 3rd ed. (Academic Press, 
2008). 

[2] Y. R. Shen, The Principles of Nonlinear Optics (John 

Wiley & Sons, 2003). 
[3] M. Kuzyk, Phys. Rev. Lett. 85, 1218 (2000). 
[4] M. Fujiwara, K. Yamauchi, M. Sugisaki, K. Yanagi, 

A. Gall, B. Robert, R. Cogdell, and H. Hashimoto, Phys. 

Rev. B 78, 161101 (2008). 
[5] J. M. Hales, J. Matichak, S. Barlow, S. Ohira, K. Yesu- 



das, J.-L. Bredas, J. W. Perry, and S. R. Marder, Science 
327, 1485 (2010). 
[6] T. Luu, E. Elliott, A. Slepkov, S. Eisler, R. McDonald, 
F. Hegmann, and R. Tykwinski, Organic letters 7, 51 
(2005). 

[7] J. C. May, I. Biaggio, F. Bures, and F. Diederich, Applied 

Physics Letters 90, 251106 (2007). 
[8] B. Champagne and B. Kirtman, Physical Review Letters 

95, 109401 (2005). 



8 



[9] B. Champagne and B. Kirtman, Journal of Chemical 

Physics 125, 024101 (2006). 
[10] D. Watkins and M. Kuzyk, JOSA B 29, 1661 (2012). 
[11] K. Tripathy, J. Moreno, M. Kuzyk, B. Coe, K. Clays, 

and A. Kelley, The Journal of chemical physics 121, 7932 

(2004). 

[12] J. Zhou, U. Szafruga, D. Watkins, and M. Kuzyk, Phys. 

Rev. A 76, 053831 (2007). 
[13] S. Shafei, M. Kuzyk, and M. Kuzyk, JOSA B 27, 1849 



(2010). 

[14] T. J. Atherton, J. Lesnefsky, G. A. Wiggers, and R. G. 

Petschek, Journal of the Optical Society of America B 

pp. 1-8 (2012). 
[15] M. Kuzyk and C. Dirk, Phys. Rev. A 41, 5098 (1990). 
[16] G. Wiggers and R. Petschek, Optics letters 32, 942 

(2007). 



